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The expansion of a classical Hamilton formalism consisting in adaptation of it to describe the 
nonequilibrium systems is offered. Expansion is obtained by construction of formalism on the ba- 
sis of the dynamics equation of the equilibrium subsystems by which the nonequilibrium system 
is represented. It has allowed removing restrictions on dynamics of the subsystems, which dic- 
tated by the requirement of monogenic and potentiality of the forces between subsystems. Modified 
Lagrange, Hamilton and Liouville equations are obtained. Some features of dynamics of nonequi- 
librium systems are considered. Connection between the equation of interaction of subsystems and 
a thermodynamic principle of energy is analyzed. 
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I. INTRODUCTION 

The attempts of further development of methods of re- 
searches of nonequilibrium systems collides with a prob- 
lem of irreversibility [1-4]. The analysis of approaches to 
solution of this problem leads us to conclusion, that the 
one of the general reason of its difficulties is a restriction 
on the frameworks of formalisms of the classical mechan- 
ics. Really, creation of the organized structures is caused 
by the dissipations. But dissipation is not present in 
Hamilton systems because the obtaining of the Hamilton 
principle based on condition that the forces of interaction 
of systems are potentially and monogenic [5]. 

Necessity of expansion of formalisms of the classical 
mechanics have arisen as a result of investigating of dy- 
namics of hard disks [6, 7]. Forces of interaction of hard 
disks systems are non-potential. Therefore it was nec- 
essary so to modify canonical Lagrange, Hamilton and 
Liouville equations that they have been applicable for 
the description of dynamics of non-potentially of inter- 
acting systems. Such modification consisted in obtain- 
ing of these equations basing on the D'Alambert prin- 
ciple under condition of non-potentiality of the forces 
[5]. From modified Liouville equation followed that non- 
potentiality of collective forces of interacting of systems 
is a necessary condition of existence of irreversible dy- 
namics [7]. 

In the nature all fundamental forces are potential. 
Therefore it was necessary to generalize results of re- 
searches of hard disks on systems of potentially inter- 
acting elements. For this purpose as a method of the de- 
scription of dynamics of nonequilibrium systems we have 
taken advantage of an opportunity of representation of 
nonequilibrium systems as a set of equilibrium subsys- 
tems (ESS). This method of system's splitting on the 
ESS is not new. It was used with success at construc- 
tion of statistical physics of equilibrium systems when 
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interaction ESS can be neglected. But it follows from a 
condition of maximum entropy that ESS for nonequilib- 
rium system will be in motion in relative to each other 
[8]. Therefore to determine their dynamics it is neces- 
sary to take into account the interaction between them. 
The equation of interaction of systems (UVS) for cal- 
culation of the ESS dynamics was obtained within the 
framework of classical mechanics [9]. Using UVS and 
basing on D'Alambert principle about equality to zero of 
a variation of work of effective forces, the modified La- 
grange, Hamilton and Liouville equations for interaction 
ESS have obtained. These equations determine expan- 
sion of a formalism of the classical mechanics, allowing 
study of the open and nonequilibrium systems. 

This paper is constructed as follows. The UVS and 
modified Lagrange, Hamilton and Liouville equations are 
obtained. The analysis of the important laws of dynamics 
of nonequilibrium systems and relation between classical 
mechanics and thermodynamics are submitted. 



II. THE EQUATION OF ONE SYSTEM 
DYNAMICS 

Let us show, how obtaining of the equations of dy- 
namics of elements and their systems are possible. The 
equation of dynamics of two interacting systems will be 
obtained in a similar way. 

Let us consider system of potentially interacting ele- 
ments. If the system is conservative, we will have for 
energy E: E = 0. The energy of the system depends 
on elements velocities, and their coordinates. Therefore 
we have: E — E{r, v), where r, v are a set of coordinates 
and velocities of the particles. In this case the equality, 
E = 0, takes place only when energy depends on two 
additive parts. One part should be a function of veloc- 
ity, and the second one - coordinates. Then energy of 
the system can be presented as E = ip [T + U] = const, 

N 

where T = ^ Ti(vi 2 ), Vi is a velocity of i- element, T is 
the kinetic energy of the system, and U(ri) is the poten- 
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tial one [10]. The (/^-function should be linear in order 
to be constant when the coordinates and velocities are 
being changed. It is always possible to represent such 
function as E = T + U by means of scale transformation 
and usage of the necessary system of coordinates. So, the 
sum of kinetic and potential energies of the system in a 
non-homogeneous space should be constant. 

Let us take a moving elementary particle with mass m 
and velocity v. The kinetic energy corresponding to the 
particle will be T(v 2 ) = mv 2 /2, and the potential energy 
-U(r), so E = mv 2 /2 + U(r) = const. In this case from 
equality E = follows that: 



v(mv + dU I dr) = 



(1) 



The eq. (1) is a balance equation of the kinetic and po- 
tential energies. As cq. (1) should be fulfilled in case of 
any direction of the vector velocity, the following require- 
ment takes place: 



mv = —dll/dr 



(2) 



It is Newton equation (NE). This equation is determining 
the connection of acceleration of the particle with the 
external force. The right hand side of eq. (2) is the active 
force. The left hand side - is the inertial force [5]. The 
particle moves along the gradient of a potential function. 
During the motion of the particle along the closed line 
in a potential field, the work of forces is equal to zero. 
Therefore dynamics of the particle is reversible. 

Let us take a system, which consists of N potentially 
interacting elements located in the non-homogeneous 
space; the mass of each element is equal to 1. The force 
acting on each element is equal to the sum of forces of 
all elements and the force caused by non-homogeneous 
space. The force between two elements is central and it 
depends on the distance between them. 

Let us represent the energy of the system as a sum 

N 

of kinetic energies of elements - TV = E mw i 2 /2, po- 

i=l 

tcntial energies in the field of external forces - UN env , 
and the potential energy of their interaction J7jv(^ij) = 

JV-1 N 

E E Uij(rij), where - is a velocity of i -element; 

i=l 

Tij = Ti — Tj - is a distance between elements i and 
j. So, E = E N + U env = T N + U N + U env = const. 
The time derivative of the energy will be as follows: 

AT N 

E = E N + U env , where E N = 5>i("M>i + £ F^); U env = 

i=l j=£i 

N 

£ ViF env. F .. = du N /dr ih Fi nv ( n ) = du env /dn. So 

i=l 

we have: 



N 

E 



ViR = 



(3) 



Where F { = mii + Fi 3 + Ff nv is effective force 
for i particle. Then the eq. (3) can be rewritten as: 



E = J2 v iFi — 0- This equality can be treated as or- 

i=l 

thogonality of the vector of effective forces with respect 
to the vector of velocities of elements of the system. If 
there are no restrictions imposed on the Vi directions, the 
requirement Fi — is satisfied [5]. Then from eq. (3) we 
obtain: 



N 



mv. 



= -J2v l F l -F? 



(4) 



i=l 



It is NE for the system's elements in non-homogeneous 
space. The eq. (4) has coincided with the corresponding 
equation in [11] though it is obtained by other way. 

Let us obtain an equation of motion of the system as 
a whole in an external field. It is known [10], that a 
motion of the system in the homogeneous space is char- 
acterized by two integrals: energy of motion of the center 
of mass and the total energy of motion of the particles 
inside of the system relative to the center of mass. We 
will name second energy as the internal energy of the 
system. Therefore it is reasonable to rewrite the eq. (3) 
having presented systems energy as the total of these 
two types of energies. For this purpose we can use the 

N 

following equality: T N = £ mv, 2 /2 = m/(2N){V^ + 



i=l 



N 



N-l N 

E E 

i=l j=i+l 

velocities of the center of mass; Rn - are coordinates 



v ij} ( a )i where Vat = — 1/N h _are 

i=l 



of the center of mass; Vi 



We will write the en- 



ergy of the system in such a way: E^ — TjJ + E 



N 

follows: 



Un- Then the eq. (3) can be written as 



T'ltr i 7712ns 
TV ' l\ 



N 

E 



v,F? 



(5) 



Where T$ = M N V N V N ;M N = mN; E^ s = T^ s + 

JV-l N 

UW S = E E Vijimvij/N + Fn). 
1=1 j=i+\ 

The term T^J* in cq. (5) determines the change of the 
kinetic energy of motion of the system as a whole, and 
the term E^ s represents change of the internal energy. 
The right hand side determines the change of energy of 
the system result of action of external forces. 

Let us to transfer to the generalized variables in eq. 
(5). For this purpose we will represent velocity of the 
motion of elements of the system as the sum of veloci- 
ties of their motion with respect to the center of mass 
of the system - ii , and velocity of the center of mass 
itself - Vn , i.e. = Vi + Vjv. Using these vari- 

N 

ables we will have the following: T/v = E , mv 2 /2 = 

1=1 

N N N 

m/{2N)V% + mV N E^+E mvf/2. As"£vi = 0, then 



3 



N N 

T N = m/(2N)V% + J2 mvf/2. Therefore £ mi 2 j2 = 

i=l i=l 

N-l N 

1/(2N) ^2 E v fj- Thus the kinetic energy of the rel- 

i=l j=i+l 

ative motion of particles of the system equals the sum of 
kinetic energies of the particles' motion with respect to 
the center of mass. 

Let us take into account that = f^ = fj — fj , where 
fi , fj - are coordinates of the elements with respect to the 
center of mass of the system. So Un {r%j ) = Un {fij ) = 

JV-l N N 

U N {fi). That is J] E v ij F ij( r ij) = E ^(n). 

i—l j—i+1 i—1 
N 

where Fi = dUN/dfi — ^ dUN/drij. By means of 

generalization of corresponding equalities for the kinetic 
and potential energies of the system, we will obtain from 
eq. (5): 

AT 

V N M N V N + mvi{ii + F(f)i) = 

N 

= -V N F env - ViFr v (R, h) (6) 

i=l 

N 

Here F 6 ™ 1 ' = £ F i nv {R, h), R-is a coordinate of ccn- 

i=l 

ter of mass. 

The eq. (6) determines the balance of energy of the 
system in non-homogeneous space. The first term in the 
left hand side corresponds to change of kinetic energy of 
motion of system as the whole. The second term deter- 
mines the change of internal energy. 

Let us take into account that F env = F env (R + fj), 
and suppose that R 3> fj. Then it is possible to ex- 
pand the force F env in a series using a small parameter, 
fi/R. Keeping the terms up to first-order of infinitesimal, 
we will have: Ff nv = Ff nv \ R + (\7Fl nv )\ R f t = F$? v + 

N N 

(VFf^ v )fi. Taking into account that E ^ — E ^ = 

i=l i=l 

N 

and J2 F to v = NF fo v = F o nV i we can set from ( 6 ) : 

i=l 

AT 

V N (M N V N ) + mvi(vi + F(f)i) w 

i=l 

AT 
i=l 

The force Fq™" is potential and depends on R. The 
right hand side of eq.(7) determines the work of external 
forces. The first term is depending on the velocity of 
systems motion and coordinate of its center of mass. It 
is determine the change of kinetic energy of system as 
the whole. 

The second term is depending on coordinates of parti- 
cles and their velocities in relative to the center of mass. 



It is determine the change of internal energy of system. 
Thus, work of external forces is splitting on two essen- 
tially different parts. 

The right hand side is the first order of infinitesimal, 
because despite R ^> f j , it value are not small. 

If the external field is homogeneous then (VF env ) i0 = 
0. The variables in the eq. (7) is separating, and we 
have: 

M N V N = F% nv (8) 

mvi = -F(fi)i (9) 

The eq. (8) is an equation of the motion of the sys- 
tem as a whole in an external field, and the eq. (9) is 
an equation of the motion of elements of the system. It 
follows from these equations, that in a homogeneous ex- 
ternal field the kinetic energy of the system as a whole 
is variable only, and the internal energy is constant. The 
change of this kinetic energy does not depend on the 
forces between elements. The system is equivalent to the 
single body with the mass which is equal to the sum of 
masses of the particles. The motion of elements of the 
system does not depend on an external field and is de- 
termined by forces of interaction only. 

Thus, the system of potentially interacting elements in 
an external field can be considered, us a structured par- 
ticle. Generally the work of external forces changes both 
its internal energy, and energy of its movement. There- 
fore dynamics of such system (or the structured particle) 
is determined by the equation (6). When the hetero- 
geneity of external forces can be neglected, dynamics of 
system is described by the Newton equation. 

III. THE EQUATION OF INTERACTION OF 
TWO SUBSYSTEMS 

Let us accept that the nonequilibrium system can be 
presented by set of moving relative to each other of ESS. 
For ESS a role of external forces carry out forces from 
other ESS. Obviously, it is impossible to neglect by het- 
erogeneity of these forces in generally Therefore dynam- 
ics ESS will be described by the equation (6) in which the 
external forces is a forces from other ESS. Hence, hav- 
ing replaced in (6) the right-hand side on these forces, 
we will obtain the UVS which describe the dynamics of 
nonequilibrium system. Let us show, how it is possible 
to obtain UVS. 

Let us the system consists of two interacting equilib- 
rium subsystems - L and K. Let us all elements to be 
identical and have the same weight 1, and L to be a num- 
ber of elements in L - subsystem, K -is a number of ele- 

L 

ments in K -subsystem, i.e. L + K = N, Vl = 1/L v i 

i=i 

K 

and Vk = 1/K J2 v i ~ are subsystems' velocities with re- 

i=l 

spect to the center of mass. Let the velocity of the center 
of mass be equal to zero, i.e. LVl + KVk = 0. 
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We can represent the energy of the system as En — 
El + E K + U mi = const, where E L and E K are the 
subsystems' energies, and U mt - is the energy of their 
interaction. According to the eq. (4), the energy of each 
subsystem can be represented as El = T L r + E L ns , Ek = 
T$ + E™, where Tf = Ml / 2 , T*T = M K V K /2, 
Ml = mL, Mk — mK. E ms - is the internal energy 
of a subsystem. The E ms consists of the kinetic energy 
of motion of the elements with respect to the center of 
mass - T ms and their potential energy - U lns , i.e. E ms = 

T ms + U ms , where U™ = E £ U iLjL (n LjL ), 

IL=1 3L=iL+l 

K-l K 

E E U iKjK (r iKjK ). The energy U mt is 



TTins 
U K 



IK = 1 jK=lK + l 



K L 

determined as U mt = £ £ Uj L j K {rj L j K ). Indexes 

jk,3L,iK,iL determine belonging of the elements to cor- 
responding subsystem. In equilibrium T tr = 0. Hence, 
if the system aspirates to equilibrium, then T tr energy 
for each subsystem will be transformed into the internal 
energy. 

We have obtained the equations of dynamics of L and 
K ESS in the following way. Let us differentiate energy 
of system on time. In order to find the equation for L - 
subsystem, at the left hand side of obtained equality we 
have kept only those terms which determine the change 
of the kinetic and potential energies of interaction of el- 
ements of L - subsystem. We replaced all other terms 
in the right hand side and combined the groups of terms 
in such a way when each group represents plurality of 
the terms forming the NE for If -subsystems elements. 
Therefore these groups are equal to zero. As a result the 
right hand side of the equation will be contain only the 
terms which determine the change of the potential energy 
of interaction of the elements L - subsystem with the el- 
ements K - subsystem. The equation for K - subsystem 
we will obtain in the same way. So we will have: 

L L-1 L L 



£ E «t 

lL=ljL=lL + l 



. f 1 . . 
L3L 1 tLjL 



E v iL F t 



K 



F ■ - 



i L = l 



K 



E v jv Fff 



JK = 1 



"OK"- 3K 



IL3K I 



FUR L ,r iK ) 



E rnv iL v 

i L = l 

(b), 

and for if -subsystem: 

K K-l K 

J2 mv lK v lK + E 

lK = l iK=ljK=iK+l 

(C). 

Here F![(R K ,r iL ) = £ *i 

3k = 1 
L 

E Fi L j K - are forces between the corresponding par- 

»l=1 

ticle of one ESS and all particles of the other ESS. The 
work of these forces determines the change of energy of 

ESS. R K = (l/K) £ r lK ; R L - (1/L) £ n L . 

i K =l ih = l 

For obtaining from the equations (b, c) the equations 
of subsystems interaction, it is necessary to present their 
left-hand side in the form of two terms. One term should 
characterize the change of energy of a motion of a ESS as 



the whole, and second is a change of its internal energy. 
With this purpose let us take into account that velocities 
of the elements can be represented as the sum of their 
velocities with respect to the center of mass of the ESS 
and the velocity of the ESS with respect to the center of 
mass of the system, Vi=Vi + V. Then with a help of (a), 
we will obtain from (b, c): 

VlMlVl+J: ]T {^[^p + 

+F iLjL ]} = -^ L -V L ^ (10) 



K 

*jr=l 3K=iK + 1 

+F iKjK ]} = $ K + V K * (11) 



$7 



= E v iL F^, 



<5> K = 



Here * = £ F k 

E 

l K = l 

The eqs. (10, 11) are equations of systems interaction 
(UVS). The first terms in the left hand side of UVS is a 
change of motion energy of ESS - "T* r ", and the second 
terms is a change of their internal energy - " E ms " . 

The right hand side of UVS are consisting of two terms 
determining the work of collective forces of interaction 
ESS. The first terms, "<I>l " and determine the 

work of the interaction force on transformation of motion 
energy of ESS into their internal energy. The second 
terms include potential force, It is a central force 

which depends on distances between the centers of mass 
of ESS: * = V(R L k), where R LK = R L - Rk- The 
work of determines the change of velocities of ESS. 
As well as in the case of the motion of the system in a 
non-homogeneous external field (see the eq. (7)), when 
Rl 3> fi L , Rk 3> ?i K this work can be estimated by 
expanding of force in a Taylor series. 

It is possible to transform the equations (10, 11) 
to the form of the equations containing forces of fric- 

L 

tion. Let us take into account that £ mv il v iL = 

t£=l 

L-1 L . L 

E E m ^L3L^L3L^ U L = E F tL^ L = 

IL = 1 jL=lL+l VL=1 

L-1 L L 

E E F iLjL v iLJL , F lL = £ dU L /dr iL . And 

»i = l 3L=iL + l 3L^IL 

for K -subsystem is similar. If also introduce the fol- 
lowing designation: a L = -(E' l [ LS + $ L )/V L 2 , a K = 

-{ET - $ K )/V K , Et s = E < (mi iL + F iL ), Ef s = 

»L = 1 

K 

£ Vi K (mvi K +Fi K ) then the eqs. (10,11) can be rewrit- 

i K = l 

ten as: 



M L V L 



- a L V L 



(12) 
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M K V K = -* - a K V K (13) 

Here coefficients "a/," , "ax" are equal to the ratio be- 
tween the energy of motion and the full kinetic energy 
of correspondent ESS. Thus, the forces of interaction of 
the ESS are divided into the potential and non-potential 
parts, and a is a similar to the friction coefficient. The 
non-potential part of the force determined by coefficient 
corresponds to dissipation. It is important to note that 
here energy is not disappeared. Actually the value of en- 
ergy is constant, but energy is redistributed between the 
elements of the subsystem. The eqs. (12, 13) are equiva- 
lent to the empirical motion equations for the structured 
particles. 

Let us perform the following replacement: Vl = 
KV L k/N, V k = -LV L k/N, where V LK = R LK =V L - 
Vk- As the center of mass of the system is immovable, 
we will have from eqs. (10, 11): 

e lk + e 'l s + ET = * (14) 

Here $ = -* L + Q K , E* K = V LK {M R V LK + tf) - is 
the change of energy of the relative motion of the ESS, 
Mr = mKL/N. 

The second and third terms in the eq. (14) determine 
the change of their internal energies. The right hand side 
determines the change of interaction energy. It is obvious 
that for an equilibrium system E L r K = 0. Then the eq. 
(14) is splitting on two independent equations for ESS. 

We see that the work of the forces between ESS goes 
not only on change of their kinetic and potential energy, 
as in the case of elementary particles. It also goes on 
change of their internal energy. Therefore the force of 
interaction ESS represents the sum of two various forces. 
The first force determines motion of ESS as the whole as 
a result of mutual transformation their kinetic and po- 
tential energies. It is potential force. The second force 
determines the work on transformation of energy of inter- 
action to the internal energy of ESS. It is non-potential 
force. This force does not change the momcntums of ESS. 

Thus, it is possible to consider equilibrium subsystems 
L and K as the structured particles. The motion of ESS 
are determined not only by transformations their kinetic 
and potential energy. But it is determined by the changes 
of their internal energy olso. Therefore the force of in- 
teraction ESS represents the sum of two various forces. 
The first force determines motion of ESS as the whole as 
a result of mutual transformation their kinetic and po- 
tential energies. It is potential force. The second force 
determines the work on transformation of motion energy 
into the internal energy of ESS. It is non-potential force. 
This force does not change the momentums of ESS. 

If L = 1, K = N — 1 the equation (10) for L - subsys- 
tem will coincide with NE. The internal energy will be 
absent, and force acting on it will be potential. Therefore 
its dynamics is reversible. Therefore reversibility of dy- 
namics of separate elements does not mean yet reversibil- 
ity of dynamics of their systems. It is essential difference 



UVS from NE for elementary particles because in accord- 
ing with NE, the elementary particles return into initial 
points if the direction of velocity of their centers of mass 
will change on opposite, but ESS does not come back into 
initial points of their phase space when reversing them 
velocities. 

Generally the potentials of interaction between parti- 
cles inside ESS can be distinct from potentials of inter- 
action of particles for different ESS. 

IV. THE LAGRANGE, HAMILTON AND 
LIOUVILLE MODIFIED EQUATIONS 

There is a question as in connection with the UVS the 
canonical Lagrange, Hamilton and Liouville are modi- 
fied. The canonical equations follow from the integral 
principle of Hamilton [5, 14]. In turn the integral prin- 
ciple of Hamilton follows from differential principle of 
D'Alambcrt in case of potentiality of active forces be- 
tween ESS. D'Alambert equation is formulated on the 
basis of NE. 

In accordance with D'Alambert principle: "the work of 
the effective forces which includes the inertial and active 
forces is equal to zero for all reversible virtual displace- 
ments compatible with the restrictions given" [5]. For 
the closed non-equilibrium system that represents them- 
selves by a set of the ESS, this principle will be as follows: 
the work of all forces of interaction of the ESS is equal 
to zero for all reversible virtual displacements compatible 
with the restrictions given. 

If the change of the internal energy of ESS can be ne- 
glected, the work on their motion will be determined only 
by the potential part of interaction ESS forces. In this 
case instead UVS we can use NE and therefore come to 
the known canonical equations of a classical mechanics 
[5] . If the neglecting by the change of internal energy is 
impossible, D'Alambert equation should be written down 
on the basis UVS which take into account the distribu- 
tion of the subsystems' motion energy as a whole between 
its internal energies in the result of the work of the non- 
potential part of forces. 

Let us take a closed non-equilibrium system consisting 
of N elements which can be represented by a set of R = 
2 an ESS. The modification of formalism is carried out 
by standard way [5]. Firstly, on the basis of UVS the 
D'Alambert equation is obtained. Based on it a principle 
of Hamilton, and then Lagrange equation are obtained. 
After that the Hamilton and Liouville equations for ESS 
are deduced. 

If there are no external restrictions on system then the 
virtual displacements of the elements and ESS can be 
combined with the real ones [5]. As virtual work of all 
forces is determined by the eq. (14) we will have: 

SW LK = SEl r K + 5ET S + 6E™ S + S$ = (15) 

The eq. (15) is D'Alambert equation for a non- 
equilibrium system. The first three terms in (15) deter- 
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mine the virtual change of energies of the ESS. The two 
last terms determine the virtual work for displacement of 
the elements of one ESS in the field of forces of the ele- 
ments of the other ESS. I.e. they determine the change 
of internal energies of the ESS. The line above 5Wlk 
means that it is just a differential form which does not 
come to a scalar function variation. 

Let us transform the eq. (15) by means of multiplying 
by dt, and integrating from t = t\ to t = t<i- As a result 
we will have: 



/2 /*2 
SW LK dt = J [6E 



\ T K + &Ef s + SE% ls + 



L K 

+ S ^F lL - Yl Sf iK F iK ]dt = (16) 

where Sr iL -is a virtual displacements of II element of L - 
subsystem; 8f iK - is a virtual displacements of %k element 
of K- subsystem. 

Let us rewrite the first three terms as follows: 

tfil E M R Vl K /2 + U LK ]dt - [ £ M R V LK 5R LK ]\\; 
8 III E rnvf L /2 + U L ]dt - [ £ mv iL 5r iL ]%; 

i L =l iL=l 

SfllZ mvfj2 + U K ]dt-[f: mv iK 5r iK )l\; 

The last two items in (16) are not the full differentials 
from any scalar function. I.e. in general it is impossible 
to represent the active forces of ESS interaction in the 
form of a gradient from a forcing function. Such forces 
are named as polygenic ones. 

Let us the virtual displacements be equal to zero at 
the ends of the interval [fi,^]- Then the eq. (16) can be 
rewritten as follows: 

SS= 8W LK dt = / { y2[d/dt(d%/dv n ) - 
J 1 J 1 „=i 

-d$/dr n + F n }}Sr n dt (17) 

N 

Here S- is action, 3 = E [ mv n + U n ] - is Lagrangian, 

n=l 

r n — r n r n ■ 

All variables are independent, and that is why in order 
to achieve simplification we have used through indexing 
for all elements of the system instead of indexing in accor- 
dance with the ESS, and at the same time new indexing 
is in accordance with the one accepted in the eq. (16). 
Thus, the action for a non-equilibrium system consists of 
the terms describing dynamics of the ESS as a whole, in- 
terior dynamics of the elements of the ESS and the term 
describing interaction of the ESS. 

It follows from the requirement of independence of the 
generalized variables that 



d/dt{d%/dv n ) - d%/dr n = -F n 



(18) 



This is the very modified Lagrange equation for a non- 
equilibrium system. The right hand side of the eq. (18) 



is equal to zero only when the relative motion of the ESS 
are absent. It should take place when the system is in 
equilibrium. 

The eqs. (10, 11) were obtained by means of indepen- 
dent methods from the requirement of the system energy 
to be constant. It allows us to write the modified La- 
grange equation for each of the ESS at once. 



d/dt(dQ L /dvi L )-dQ L /dn 



-F; 



K 



(19) 



d/dt(d$S K /dv iK ) - d$ K /dr lK = -Ft K (20) 

Here 5l and 3^- - are Lagrange functions for the ESS. 
Now let's obtain modified Hamilton equation. For 3 

N 

the differential can be written as cRy= [(d^s / dv n )dv n + 

71=1 

{dQ/dr n dr n )} + d$s/dtdt. With the help of Leg- 

N 

endry transformation, we will obtain: d[ Y^, Pn,v n — 

n=l 

N 

9] = E [v n dpn ~ (d$/dr n )dr n ] - (d$s/dt)dt. Since 

n=l 

N 

d^s/dt~—dH/dt , where H = d[ ^ p n v n — we will 

n=l 

have: 



dH/dr n = -p n - F n 



dH/dp ri 



(21) 



(22) 



The above are modified Hamilton equations for a non- 
equilibrium system. The right hand side of the eq. (21) 
denotes non-potential forces. Modified Hamilton's equa- 
tions for the ESS can be obtained in the same way. 

Now let us obtain Liouville equations for the ESS and 
full system. For this porpoise let us take a generalized 
current vector Jl = JL{v n ,Pn) of the L-subsystem (here 
"n" index corresponds only to the particles of the L- 
subsystcm). From the eqs. (21, 22) we find: divJi, = 
L L 

J2 (dv n /dr n + dp n /dp n ) = - Y, dF*/dp n . The conti- 



71=1 



71=1 



nuity equation is a differential form of the particle num- 
ber conservation law: dfL/dt + div(fi J£) = 0, where 
/t = /l (r, p, t) is a normalized distribution function for 
the elements of the L- subsystem. With the help of the 
continuity equation we will have: dfL/dt = dfL/dt + 

L 

E {v n df L /dr n + p n df L /dp n )=df L /dt + div{f L J L ) - 

71=1 

L 

f L divJ L = -fLdivJ L =fL E dF^/dp n . Thus we have: 

71 = 1 



df L /dt = f L V dF^/dp n 



(23) 



71 = 1 
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The eq. (23) is modified Liouville equation for the L 
-subsystem. 

There is a formal solution of the eq. (23): /z, = 

L 

constxexpf{ E dF^ /dp n }dt. From this solution fol- 

n=l 

lows that the distribution function of a subsystem de- 
pends on time if force acting on the subsystem depends 
on velocities of the elements of the subsystem. It is sim- 
ilar result as for a hard disks systems [7] . 

The modify Liouville equation for system can be obtain 
from cqs. (21, 22) if n = 1, 2...N. We will have: 

df N /dt = f N dF n /d Pn (24) 

Despite absence of external forces for a system, the right 
hand side of the eq. (24) is not equal to zero. It is caused 
by the fact that the energy of subsystems' interaction 
does not include into ESS energy. Therefore the right 
hand side of the eq. (24) will not be equal to zero until 
the energy of interaction of the ESS is not transformed 
into their internal energy. 

The reason of distinction between two descriptions of 
system (in a framework of canonical and modified equa- 
tion) is the total work of potential forces between par- 
ticles in the closed system is equal to zero. But the 
work on moving ESS is not equal to zero. The New- 
ton equation for select particles can not determine this 
work. Therefore modified Lagrange, Hamilton and Liou- 
ville equations as against their canonical prototypes, are 
applicable for the description of dynamics of noncquilib- 
rium systems. 



V. SOME PROPERTIES OF DYNAMICS OF 
INTERACTING SUBSYSTEMS 

In accordance with UVS, in order to describe a non- 
equilibrium system, the energy should be divided into 
T tr and E ms . As UVS includes the parameter describ- 
ing the degree of non-equilibrium of a system, it allows 
studying evolution of the system. Let us to analyze some 
important cases of the system's dynamics. 

In general case the non-potential forces can be 

L 

rewritten as: $ L = E $ iL (R K ,r iL ,v iL ), § K = 

»£=1 

K 

E &i K {RL,r iK ,Vi K ). If the sizes of ESS are less 
than the distances between them, we will have &l ~ 
E Fg(R K )v iL = F* E < = 0, and <S> K „ 

i L = l il = 1 

E Ff K {R L )v iK = Ft K E v iK = 0. In this case the 

■ k = 1 iff=l 

eqs. (10, 11) are the Newton equations for two hard 
bodies. This case is similar to the case of homogeneity 
of an external field of forces. 

In the case when = 0, Vl = Vk = 0, the first terms 
of the left side of UVS are equal to zero and the system is 



in equilibrium. The full energy of the system is equal to 
the sum of E ms for each subsystem. Such system can be 
studied with the help of canonical Hamilton equations. 
ESS. Gibbs used this approach for creation of statistical 
mechanics of equilibrium systems [3, 11]. 

Now let us velocities of the elements with respect to 
the centers of mass of corresponding ESS be equal to 
zero during interaction of the ESS (this is a hard body 
approach). Then the right hand side, and also the second 
terms of left hand side of UVS are equal to zero, and UVS 
is transformed into NE for two hard bodies. 

It is possible to come from classical mechanics to ther- 
modynamics with the help of the UVS. Really, the right 
hand side of eqs. (10, 11) determines an change of ESS 
energy. The first term of the left hand side of UVS de- 
termines the change of the motion energy of ESS as the 
whole. In thermodynamics it corresponds to work, which 
is carried out by external forces acting on ESS on the part 
of an environment. The second term of the left hand side 
corresponds to increase in the entrance energy of a ESS. 
In thermodynamics this term corresponds to the change 
of thermal energy of ESS. 

Let us consider the relation between the UVS and 
the basic equation of thermodynamics in details. We 
can write the basic equation of thermodynamics [8,13]: 
dE = dQ — PdY. Here, according to common termi- 
nology, E is energy of a ESS; Q is thermal energy; P 
is pressure; Y is volume. We take into account that 
N >> l,Mjv = 1. The energy change of the selected 
ESS is due to the work made by external forces. There- 
fore, the change in full energy of a subsystem corresponds 
to dE. 

The change of kinetic energy of motion of a ESS as 
the whole, dT tr , corresponds to the term PdY. Really, 
dT tr = VdV = VVdt = Vdr = PdY 

Let us determine, what term in UVS corresponds to 
the change of the internal energy of ESS. As follows from 
virial theorem [10], if the potential energy is a homoge- 
neous function of second order of the radiuses- vectors, 
then E ins = 2T ins = 2U ms . The line denotes the time 
average. Earlier we obtained that the energy, E lns , in- 
creases due to T tr . But the opposite process is impos- 
sible. Therefore the change of the term Q in the UVS 
corresponds to the change of the energy E ms . 

Let us consider the ESS near to equilibrium. If 
the ESS consist of N elements, the average energy of 
each element becomes, E ins = E ins /N = kTq™ s . Now 
let the internal energy increases with dQ. According 
to the virial theorem, keeping the terms of the first 
order, we have: dQ w T£ ns [dE ins /Tjj ns ] = T^ ns [dv/v ], 
where vo is the average velocity of an element, and 
dv is its change. For ESS in equilibrium, we have 
dv/vo <~ dT/T, where T is the phase volume of a ESS, 
dT will increase due to increasing of the ESS energy on 
the value, dQ. By keeping the terms of the first or- 
der we get: dQ w T l ns dY/Y = T l ns d\nT. By definition 
d\nT = dS ms , where S ms is a subsystem entropy [9]. So, 
near equilibrium we have dQ w T^ ns dS ms . So formally 



8 



according with UVS the entropy production in the non- 
equilibrium system is determined by transformation of 
the motion energy of ESS into the internal energy. If it so 
eventually relative velocities of ESS go to zero. In result 
the energy of relative motion of ESS will be completely 
transforms into the internal energy and the systems equi- 
librates. It means that energy of motion of a ESS goes on 
increase of entropy. Therefore the deviation of entropy 
from equilibrium can be determined by the next formula 
[7]: 

AS = £^£/£%^} (25) 

1=1 k=l J s 

Here E rn ' is the kinetic energy of ESS; to; is the number 
elements in ESS "Z"; R is the number of ESS; s is num- 
ber of the external element which collided with internal 
element k; F™ s l is a force, acted on element fc-elemcnt; 
Vk -is a velocity of the k- element. 

The formula (25) means, that a entropy production 
is provided by the energy of relative motion of ESS. It 
corresponds to entropy definition for the non-equilibrium 
systems, which was offered in [8]. 

But for informal connection of the equations of interac- 
tion of systems with thermodynamics, the increase of in- 
ternal energy due to the energy of relative motion should 
be irreversible. Below some arguments for the benefit of 
such increase will be offered. 

In connection with the law of momentum preservation 
of ESS the internal energy of the ESS is gradually in- 
creased due to energy E l £ K . I- e - there is a gradual de- 
crease of Vlk- To explain this fact let us take an equi- 
librium system and divide it into ESS. It follows from 
the requirement of equilibrium that the forces between 
the ESS are equal to zero. But if there are no external 
forces with respect to a subsystem, the subsystem can- 
not start moving. Indeed in this case the subsystem can 
start moving only due to its internal energy. But it is 
prohibited by the law of conservation of momentum of 
the subsystem. I.e. if the system is in equilibrium then 
it will be always remain in equilibrium state. 

Now let us consider two interacting ESS which are in 
relative motion. The relative velocity of ESS can not be 
increased due to the energy of interaction of the ESS as 
this energy is itself determined by the relative velocities 
of the ESS. As well as in the previous case, this velocity 
can not be increased due to the internal energy of the 
ESS. So relative velocity can be decreased only. 

It is follows from these arguments that UVS character- 
izes motion of system to an equilibrium state as a result 
of transition of energy of relative motion of ESS into in- 
ternal energy. When system is in equilibrium, or when it 
is possible to freeze internal degrees of freedom of ESS, 
i.e. to neglect the change of their internal energy, then 
UVS is reduced to NE. Hence, it is possible to assert, 
that the broken of the time-symmetry of the systems dy- 
namics is caused by the increasing of the internal energy. 
The decreasing of internal energy is prohibited by the law 



of conservation of momentum of the subsystem. 



VI. CONCLUSION 

Expansion of a formalism of classical mechanics and 
approach to the analysis of nonequilibrium systems are 
based on the key statement that closed non-equilibrium 
systems can be represented as a set of moving ESS [8]. 
This statement allowed transforming the tasks of dynam- 
ics of a non-equilibrium system to the tasks of dynam- 
ics of ESS. The offered approach has allowed removed 
the restrictions of applicability formalism of the classical 
mechanics, bound with the requirement of conservatism 
of systems. It has appeared enough for elimination of 
contradictions between a classical mechanics and ther- 
modynamics. 

Having expressed energy of system in the form of en- 
ergy of a motion of ESS as the whole, their internal ener- 
gies and energy of their interaction and differentiated it 
on the time, the UVS has been obtained. This equation 
is equivalent to the equation of structured particles. 

In agreement with UVS the dynamics of ESS is deter- 
mined by the collective force of their interaction. This 
force consists of two parts. The potential part is deter- 
mining the transformation of the interaction energy into 
the motion energy of ESS. Non-potential force is deter- 
mining the transformation of the ESS motion energy into 
internal energy. Thus, unlike dynamics of an elementary 
particle, the dynamics of ESS is determined by the work 
of interaction forces which will transform energy of inter- 
action not only into potential and kinetic energy of ESS 
motion as the whole, how we have in the case of elemen- 
tary particle, but also into the internal energy. Therefore 
based on NE which is valid only for elementary particles, 
it is impossible to describe the dynamics of the nonequi- 
librium system. 

Based on UVS the modified Lagrange, Hamilton and 
Liouville equations was obtained. These equations arc 
applicable to the study of dynamics of the open non- 
equilibrium systems because its have allowed taking into 
account the work of non-potential part of force on trans- 
formation of energy of ESS motion into their internal 
energy. 

The following mechanism of irreversibility can be pro- 
posed: energy of ESS relative motion is transformed into 
their internal energy as a result of the work of the non- 
potential part of force of interaction of ESS. The decreas- 
ing of internal energy is prohibited by the law of conser- 
vation of momentum of the ESS. Therefore the system 
equilibrates, when relative motion of ESS is absent. 

The explanation of the First law of thermodynamics is 
based on the fact that the work of subsystems' interaction 
forces changes both the energy of their motion and their 
internal energy. The explanation of the Second law of 
thermodynamics is based on the condition of irreversible 
transformation of the subsystems' relative motion energy 
into their internal energy. 
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The solution many problems are necessary for the fur- 
ther development of the offered approach to studying 
nonequilibrium systems and the strict proof of the con- 
clusions offered here. Questions about the ways of split- 
ting of nonequilibrium system on ESS, about character 



of transformation of energy of relative motion of ESS in 
their internal energy, etc., are the most important be- 
tween them. 
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